x’ frontiers | Frontiers in Astronomy and Space Sciences 


ORIGINAL RESEARCH 
published: 08 July 2022 
doi: 10.3389/fspas.2022.876870 


OPEN ACCESS 


Edited by: 

André Canosa, 

UMR6251 Institut de Physique de 
Rennes (IPR), France 


Reviewed by: 

Audrey Coutens, 

UMR5277 Institut de recherche en 
astrophysique et planétologie (IRAP), 
France 

Donghui Quan, 

Eastern Kentucky University, 

United States 


*Correspondence: 
Victor M. Rivilla 
vrivilla@cab.inta-csic.es 


Specialty section: 

This article was submitted to 
Astrochemistry, 

a section of the journal 

Frontiers in Astronomy and Space 
Sciences 


Received: 15 February 2022 
Accepted: 02 June 2022 
Published: 08 July 2022 


Citation: 

Rivilla VM, Jiménez-Serra |, 
Martin-Pintado J, Colzi L, Tercero B, 
de Vicente P, Zeng S, Martin S, 
Garcia de la Concepcién J, 
Bizzocchi L, Melosso M, Rico-Villas F 
and Requena-Torres MA (2022) 
Molecular Precursors of the RNA- 
World in Space: New Nitriles in the 
G+0.693-0.027 Molecular Cloud. 
Front. Astron. Space Sci. 9:876870. 
doi: 10.3389/fspas.2022.876870 


Check for 
updates 


Molecular Precursors of the 
RNA-World in Space: New Nitriles in 
the G+0.693-0.027 Molecular Cloud 


Victor M. Rivilla'*, Izaskun Jiménez-Serra', Jesus Martin-Pintado', Laura Colzi', 

Belén Tercero”, Pablo de Vicente*, Shaoshan Zeng*, Sergio Martin*®, 

Juan Garcia de la Concepcion", Luca Bizzocchi’, Mattia Melosso”®, Fernando Rico-Villas' 
and Miguel A. Requena-Torres*'° 


'Centro de Astrobiologia (CSIC-INTA), Madrid, Spain, 7Observatorio Astronémico Nacional (OAN-IGN), Madrid, Spain, 
°Observatorio de Yebes (OY-IGN), Guadalajara, Spain, “Star and Planet Formation Laboratory, Cluster for Pioneering Research, 
RIKEN, Saitama, Japan, °European Southern Observatory, ALMA Department of Science, Santiago, Chile, SJoint ALMA 
Observatory, Department of Science Operations, Santiago, Chile, ‘Department of Chemistry “Giacomo Ciamician”, University of 
Bologna, Bologna, Italy, "Scuola Superiore Meridionale, Universita di Napoli Federico II, Naples, Italy, (Department of Astronomy, 
University of Maryland, College Park, ND, United States, "Department of Physics, Astronomy and Geosciences, Towson 
University, Towson, MD, United States 


Nitriles play a key role as molecular precursors in prebiotic experiments based on the RNA- 
world scenario for the origin of life. These chemical compounds could have been partially 
delivered to the young Earth from extraterrestrial objects, stressing the importance of 
establishing the reservoir of nitriles in the interstellar medium. We report here the detection 
towards the molecular cloud G+0.693-0.027 of several nitriles, including cyanic acid 
(HOCN), and three C4H3sN isomers (cyanoallene, CHsCCHCN; propargyl cyanide, 
HCCCHsCN; and cyanopropyne (CH3CCCN), and the tentative detections of 
cyanoformaldehyde (HCOCN), and glycolonitrile (HOCH2CN). We have also performed 
the first interstellar search of cyanoacetaldehyde (HCOCH>2CN), which was not detected. 
Based on the derived molecular abundances of the different nitriles in G+0.693-0.027 and 
other interstellar sources, we have discussed their formation mechanisms in the ISM. We 
propose that the observed HOCN abundance in G+0.693-0.027 is mainly due to surface 
chemistry and subsequent shock-induced desorption, while HCOCN might be mainly 
formed through gas-phase chemistry. In the case of HOCH2CN, several grain-surface 
routes from abundant precursors could produce it. The derived abundances of the three 
C4H3N isomers in G+0.693-0.027 are very similar, and also similar to those previously 
reported in the dark cold cloud TMC-1. This suggests that the three isomers are likely 
formed through gas-phase chemistry from common precursors, possibly unsaturated 
hydrocarbons (CH3CCH and CHsCCH>) that react with the cyanide radical (CN). The rich 
nitrile feedstock found towards G+0.693-0.027 confirms that interstellar chemistry is able 
to synthesize in space molecular species that could drive the prebiotic chemistry of the 
RNA-world. 
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1 INTRODUCTION 


Life on Earth appeared about 3.8 billion years ago, around 700 
Myr after the formation of the planet (Pearce et al., 2018), but we 
still do not know the mechanisms that made it possible. One of 
the most supported hypotheses for the origin of life is known as 
the RNA world (Gilbert 1986), in which RNA could have 
performed both metabolic and genetic roles. The process by 
which inert matter generated first the building blocks of RNA, 
ribonucleotides, and ultimately RNA itself, remains a mystery. 
Recent laboratory experiments mimicking prebiotic conditions 
have shown that ribonucleotides could be synthesized starting 
from simple molecules (e.g. Powner et al., 2009; Patel et al., 2015; 
Becker et al., 2019). A plausible origin of this prebiotic material is 
extraterrestrial delivery (Ord 1961; Chyba and Sagan 1992; 
Cooper et al., 2001) during the heavy bombardment of 
meteorites and comets that occurred around 3.9 billions ago 
(Marchi et al., 2014). These basic molecular precursors may have 
been already formed prior to the formation of the Solar System, in 
its parental molecular cloud, through the chemistry that takes 
place in the interstellar medium (ISM). Therefore, the study of the 
molecular complexity of the ISM can provide us an illustrative 
view of the chemical reservoir that could have contributed to feed 
the prebiotic chemistry on the primitive Earth, and could 
potentially develop similar processes in other places in the 
Galaxy under favourable Earth-like planetary environments. 

In the last decades, and especially in the last years, 
astrochemistry has shown that interstellar chemistry is able to 
synthesize building blocks of key biomolecules. Several of the 
precursors of ribonucleotides spotted by the prebiotic 
experiments in the laboratory have been detected in the ISM, 
like cyanoacetylene (HC3N, Turner 1971), cyanamide (NH2CN, 
Turner et al., 1975), glycolaldehyde (CH,OHCHO, Hollis et al. 
2004, urea (NH,CONH), Belloche et al. 2019), hydroxylamine 
(NH,OH, Rivilla et al. 2020), and 1,2-ethenediol ((CHOH),; 
Rivilla et al. 2022a). Among the key simple molecular 
precursors required for the RNA world, numerous works 
have stressed the dominant role of a particular family of 
compounds, nitriles, which are molecules with the C*N 
moiety. This simple but highly versatile functional group 
offers a unique potential to build-up molecular complexity 
and activate efficiently the formation of ribonucleotides 
(Powner et al. 2009; Powner and Sutherland 2010; Patel 
et al. 2015; Mariani et al. 2018; Becker et al. 2019; Menor 
Salvan et al. 2020), and also other key biomolecules such as 
peptides or nucleobases (Menor-Salvan and Marin-Yaseli 2012; 
Canavelli et al. 2019; Foden et al. 2020). 

With the aim of extending our knowledge on the chemistry of 
nitriles in the ISM, in this work we have searched for more nitriles 
towards the molecular cloud G+0.693-0.027 (hereafter G+0.693), 
including some with increasing complexity that have been 
proposed as important precursors of prebiotic chemistry. This 
cloud, located in the Sgr B2 region of the center of our Galaxy, is 
one of the most chemically rich sources in the ISM. Numerous 
nitrogen-bearing species, including nitriles, have been detected 
(see Zeng et al., 2018; Rivilla et al., 2019b, 2021b): cyanoacetylene 
(HC3N), acetonitrile (CH3CN), cyanamide (NH,CN), the 


cyanomethyl radical (H,CCN), cyanomethanimine 
(HNCHCN), and the cyanomidyl radical (HNCN). In this 
work we report the detection of cyanic acid (HOCN), the 
tentative detections of glycolonitrile (HOCH,CN) and 
cyanoformaldehyde (HCOCN), and the first interstellar search 
of cyanoacetaldeyde (HCOCH2CN) in the ISM, for which we 
provide an abundance upper limit. We have also searched for 
three unsaturated carbon-chain nitriles, the CjH3N isomers. We 
report the detection of cyanopropyne (CH3;CCCN), and the 
second detections in the ISM of cyanoallene (CH,CCHCN) 
and propargyl cyanide (HCCCH,CN), detected previously only 
towards the TMC-1 dark cloud (Lovas et al. 2006; McGuire et al. 
2020; Marcelino et al. 2021). In Section 2 we present the data of 
the observational survey, in Section 3 we describe the line 
identification and analysis, and present the results of the line 
fitting, and in Section 4 we discuss about the interstellar 
chemistry of the different species and their possible roles in 
prebiotic chemistry. 


2 OBSERVATIONS 


A high sensitivity spectral survey was carried out towards G+0.693. 
We used both IRAM 30 m telescope (Granada, Spain) and Yebes 
40 m telescope (Guadalajara, Spain). The observations were centred 
at o(J2000.0) = 1774722', and 6(J2000.0) = — 28°21'27". The 
position switching mode was used in all the observations with 
the off position located at Aw = — 885”, Ad = 290” from the 
source position. During the IRAM 30m observations the dual 
polarisation receiver EMIR was connected to the fast Fourier 
transform spectrometers (FFTS), which provided a channel width 
of 200 kHz. In this work we have used data covering the spectral 
windows from 71.8 to 116.7 GHz, 124.8-175.5 GHz, and 199.8 — 
238.3 GHz. The spectra were smoothed to velocity resolutions of 
1.0 — 2.6 km s ', depending on the frequency. The observations with 
the Yebes 40m radiotelescope used the Nanocosmos Q-band 
(7mm) HEMT receiver (Tercero et al., 2021). The receiver was 
connected to 16 FFTS providing a channel width of 38 kHz and a 
bandwidth of 18.5 GHz per polarisation, covering the frequency 
range between 31.3 and 50.6 GHz. The spectra were smoothed to a 
resolution of 251 kHz, corresponding to velocity resolutions of 1.5 — 
2.4km s '. The noise of the spectra depends on the frequency range, 
with values in antenna temperature (T4 ) as low as 1.0 mK, while in 
some intervals it increases up to 4.0 — 5.0 mK, for the Yebes data, and 
1.3 to 2.8 mK (71 — 90 GHz), 1.5 to 5.8 mK (90 — 115 GHz), ~10 mK 
(115 — 116 GHz), 3.1 to 6.8 mK (124 — 175 GHz), and 4.5 to 10.6 mK 
(199 — 238 GHz), for the IRAM 30m data. The line intensity of the 
spectra was measured in units of TX as the molecular emission 
toward G+0.693 is extended over the beam (Requena-Torres et al., 
2006, 2008; Zeng et al., 2018, 2020). 


3 ANALYSIS AND RESULTS 


Figure 1 shows the nitriles analysed in this work, which include four 
oxygen-bearing nitriles: cyanic acid (HOCN), cyanoformaldehyde 
(or formyl cyanide, HCOCN),  glycolonitrile (or  2- 
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FIGURE 1 | Three-dimensional representation of the oxygen-bearing nitriles (upper panel) and the three C4HsN isomers (lower panel) analysed in this work. White, 
gray, red, and blue corresponds to hydrogen, carbon, oxygen and nitrogen atoms, respectively. 


HOCH2CN 
glycolonitrile 


HCOCH2CN 
cyanoacetaldehyde 
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TABLE 1 | Spectroscopy of the molecules analysed in this work. 


Molecule Catalog Entry Date Line List Dipole Moment 
ref ref 

HOCN CDMS 43510 May 2009 Brunken et al. (2009) Bruinken et al. (2009) 

HCOCN CDMS 55501 June 2006 Bogey et al. (1995) Csaszar (1989) 

HOCH2CN CDMS 57512 March 2017 Margulés et al. (2017) Margulés et al. (2017) 

HCOCH2CN MADCUBA - January 2022 Mellendal et al. (2012) Mollendal et al. (2012) 

HCCCH2CN CDMS 65514 September 2020 Demaison et al. (1985) Demaison et al. (1985) 
McGuire et al. (2020) 

CHsCCHCN CDMS 65506 March 2006 Bouchy et al. (1973) Bouchy et al. (1973) 
Schwahn et al. (1986) 

CH3CCCN CDMS 65503 April 2004 Moises et al. (1982) Bester et al. (1984) 


Bester et al. (1983) 


The molecular catalog, number and date of the entry, and the references for the line lists and dipole moments are listed. 


hydroxyacetonitrile, HCOCH,CN), cyanoacetaldehyde (or 3- 
oxopropanenitrile, HCOCH,CN); and three C,H3;N isomers: 
cyanoallene (or 2,3-butadienenitrile, CH,CCHCN), propargyl 
cyanide (or 3-butynenitrile, HCCCH,CN), and cyanopropyne (or 
2-butynenitrile, CH3;CCCN). The identification and fitting of the 
molecular lines were performed using the SLIM (Spectral Line 
Identification and Modeling) tool within the MADCUBA 
package’ (version 09/11/2021; Martin et al, 2019). SLIM 
generates synthetic spectra under the assumption of Local 
Thermodynamic Equilibrium (LTE), using the spectroscopy 
provided by laboratory experiments assisted by theoretical 
calculations. Table 1 lists the spectroscopic references of all the 
molecules analysed. We have used entries from the Cologne 


‘Madrid Data Cube Analysis on ImageJ is a software developed at the Center of 
Astrobiology (CAB) in Madrid; https://cab.inta-csic.es/madcuba/. 


Database for Molecular Spectroscopy (CDMS, Endres et al. 
2016), which are based on the laboratory works and theoretical 
calculations indicated in Table 1. Moreover, we implemented into 
MADCUBA the spectroscopy of HCOCH,CN from Mollendal et al. 
(2012). 

To evaluate if the molecular transitions of interest are blended 
with emission from other species, we have also considered the 
LTE model that includes the total contribution of all the 
species that have been identified so far towards G+0.693 
(e.g., Requena-Torres et al., 2008; Zeng et al., 2018; Rivilla 
et al., 2019a, 2020; Jiménez-Serra et al., 2020; Rivilla et al., 
2021a,b; Zeng et al., 2021; Rodriguez-Almeida et al., 2021a,b; 
Rivilla et al., 2022a,b). To derive the physical parameters of 
the molecular emission, we used the AUTOFIT tool of SLIM, 
which finds the best agreement between the observed spectra 
and the predicted LTE model, and provides the best solution 
for the parameters, and their associated uncertainties (see 
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FIGURE 2| Selected cyanic acid (HOCN) transitions (see Table 2) detected towards the G+0.693 molecular cloud. The best LTE fit derived with MADCUBA for the 
HOCN emission is shown with a red curve, while the blue curve shows the total emission considering all the species identified towards this molecular cloud. The y-axis 
shows the line intensity in antenna temperature scale (T4) in Kelvin, and the x-axis shows the frequency in GHz. 


details of the formalism used in Martin et al., 2019). The free 
parameters of the model are: molecular column density (N), 
excitation temperature (T.x), linewidth (or full width at half 
maximum, FWHM), and velocity (visr). We have left these 
four parameters free whenever possible, providing their 
associated uncertainties. For the cases in which the 
algorithm used by AUTOFIT does not converge, we have 
fixed some of them to allow the algorithm to converge. In the 
following, we present the analysis of the different molecules 
studied. For each species, we have applied AUTOFIT using 
unblended transitions and transitions that, while partially 
blended with other species already identified in G+0.693, 
properly reproduces the observed spectra. We note that for 
all molecules the transitions that are not shown are consistent 
with the observed spectra, but they are heavily blended with 
lines from other molecular species or they are too weak to be 


detected, according to the line intensities predicted by the 
LTE model. 


3.1 Oxygen-Bearing Nitriles 


3.1.1 Cyanic Acid (HOCN) and Cyanoformaldehyde 
(HCOCN) 


HOCN was already reported towards G+0.693 by Briinken et al. 
(2010) (their source Sgr B2 (20,100)”), and also by Zeng et al. 
(2018) using in both cases less sensitive observations. We provide 
here a new analysis using deeper observations. We have detected 
six transitions of this species that are completely unblended, which 


?The position of this source is offset in (a, 5) by (20”, 100”) with respect to that of 
Sgr B2(M), see Briinken et al. (2010). 
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TABLE 2 | List of detected transitions of the oxygen-bearing nitriles analysed in this work. We indicate the frequency, quantum numbers, logarithm of the Einstein coefficient 
(Au), energy of the upper levels of each transition (E,), and information about the possible blending by other identified or unidentified (U) species towards G +0.693. 


Molecule Frequency (GHz) Transition Jx,K, 
HOCN 41.9508371 20,2 - 10,1 
HOCN 83.5383960 41.4 - 31,3 
HOCN 83.9005702* 4o,4 — 30,3 
HOCN 84.2524547 Aye— Si 
HOCN 104.8746777° 50,5 - 404 
HOCN 125.8480951 60,6 — 50,5 
HOCN 146.82068467 70,7 - 606 
HOCN 167.7923140° 80,8 - 70,7 
HCOCN 72.1192555 14.4 - 00,0 
HCOCN 81.433113 24,2 - 101 
HCOCN 90.5710141 314.3 — 202 
HCOCN 99.5348108 44.4 - 30,3 
HCOCN 108.3274964 54.5 — 4o,4 
HCOCN 207.3132961 220 - 14,1 
HOCH2CN 35.934379 44.4 - 31,3 
HOCH2CN 37.781601 41.3 - 310 
HOCH2CN 44.907438 5is—4ia 
HOCH,CN 45.975528 505 — 40.4 
HOCH2CN 75463333 817-716 


“Transition detected in Zeng et al. (2078). 


Transition tentatively detected in Zeng et al. (2018). 


logA.y (s~*) 


5.3239 
—4.4087 
—4.3750 
—4.3976 
—4.0746 
—3.8304 
3.6248 
3.4472 


—5.0767 
—4.9642 
—4.8461 
—4.7342 
—4.6301 
3.7478 


—5.8654 
—5.8001 
5.5549 
-5.5069 
—4,8528 


E, (K) 


3.0 
42.2 
10.1 
42.3 
15.1 
21.4 
28.1 
36.2 


3.5 
44 
5.7 
7.6 
9.8 
13.4 


5.7 
5.9 
7.9 
6.6 
17.7 


Blending 


unblended 

blended with U 
unblended 
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u 
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FIGURE 3 | Selected cyanoformaldehyde (HCOCN) transitions (see Table 2) detected towards the G+0.693 molecular cloud. The best LTE fit derived with 
MADCUBA for the HCOCN emission is shown with ared curve, while the blue curve shows the total emission considering all the species identified towards this molecular 
cloud. The y-axis shows the line intensity in antenna temperature scale (T,') in Kelvin, and the x-axis shows the frequency in GHz. 
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FIGURE 4| Selected transitions of glycolonitrile (HOCH2CN; see Table 2) detected towards the G+0.693 molecular cloud. The best LTE fit derived with MADCUBA 


for the HOCH2CN emission is shown with a red curve, while the blue curve shows the total emission considering all the species identified towards this molecular cloud. 
The y-axis shows the line intensity in antenna temperature scale (74) in Kelvin, and the x-axis shows the frequency in GHz. 


75.44 75.46 


are shown in Figure 2, and listed in Table 2. These transitions 
include the three transitions identified by Zeng et al. (2018), the 
confirmation of the 89g — 7o,7 transition tentatively detected in that 
work (see their Figure B15), and two new transitions (Table 2). The 
best LTE fit derived by MADCUBA, where all parameters were left 
free, is shown in Figure 2, and the derived physical parameters are 
presented in Table 4. We obtained a column density of 
(2.13+ 0.04) x 10!3 cm (Table 4), which translates into a 
molecular abundance with respect to molecular hydrogen of 
1.6 x 10°'°, using N(H2) = 1.35 x 10°? cm™ from Martin et al. 
(2008). The results are consistent, within the uncertainties, with 
those derived by Zeng et al. (2018). 

We also report here the first tentative detection of HCOCN 
towards G+0.693. Figure 3 shows that the 3,3; - 292 
(90.5710141 GHz) and 4,4 - 393 (99.5348108 GHz) transitions 
are unblended, while other transitions are partially blended with 
other species (Table 2). To perform the fit, we fixed T.,, FWHM, and 
Visp to the ones derived from HOCN. We obtained a column density 
of (0.76+ 0.11) x 10! cm”, almost one order of magnitude lower 
than the upper limit reported by Zeng et al. (2018) of <6 x 10% 
cm™ towards G+0.693. The derived molecular abundance is 6 x 
10’, which is very similar to that found in the TMC-1 dark cloud 
by Cernicharo et al. (2021). The HOCN/HCOCN ratio is ~2.8. 


3.1.2 Glycolonitrile (HOCHsCN) 


This species is also tentatively detected towards G+0.693. We 
show in Figure 4 two molecular transitions of HOCH,CN that 


are unblended (Table 2), and those partially blended with other 
species already identified in this cloud. To perform the fit, we 
fixed T,, and FWHM to the ones derived for HOCN, and used 
Visr = 67 km s"', which best reproduces the velocity of the two 
unblended transitions. We obtained a column density of 
(0.8+ 0.2) x 10!3 cm™? (Table 4), and a molecular abundance 
of 6 x 10 "’, very similar to that of HOCN. 


3.1.3 Cyanoacetaldehyde (HCOCH2CN) 


This molecule is not currently included in any of the commonly 
used molecular databases such as CDMS or the Jet Propulsion 
Laboratory catalog (JPL; Pickett et al., 1998). The conformational 
energy landscape of HCOCH2CN and the effects of the large 
amplitude motions on its rotational spectrum have been 
described in detail by Mollendal et al. (2012). We have used 
the spectroscopic information provided in this work to 
implement it into MADCUBA. The most stable rotamer 
(referred to as species I in the cited reference) possesses two 
equivalent positions in the electronic energy potential function 
for rotation about its C,;-C, bond (see Figure 1 of Mollendal 
et al., 2012). They are separated by a barrier of 0.84kJ mol ' 
(computed at MP2 level) at the exact antiperiplanar 
conformation. Large amplitude vibrations and tunneling for 
the torsion about the C,-C, bond leads to the existence of 
two closely spaced energy levels for the ground state labelled 
with a plus sign (+) for the lowest-energy level and with a minus 
sign (—) for the higher-energy level. These two states are separated 
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FIGURE 5 | (continued). 


by an energy difference AE/h of ~58.8 GHz. For the present reported by Mollendal et al. (2012) using the same set of 
spectral calculation we have reanalysed the rotational transitions spectroscopic parameters employed in their fit 1 (see their 
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FIGURE 5 | Selected transitions of cyanoallene (CH2CCHCN; see Table 3) detected towards the G+0.693 molecular cloud. The best LTE fit derived with 
MADCUBA for the CH2CCHCN emission is shown with a red curve, while the blue curve shows the total emission considering all the species identified towards this 
molecular cloud. The y-axis shows the line intensity in antenna temperature scale (T,‘) in Kelvin, and the x-axis shows the frequency in GHz. 
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TABLE 3 | List of observed transitions of the C,H3N isomers analysed in this work. We indicate the frequency, quantum numbers, Einstein coefficient (A.), energy of the 
upper levels of each transition (E,), and information about the possible blending by other identified or unidentified (U) species towards G +0.693. 


Molecule 


HsCCHCN 
H2CCHCN 
HoCCHCN 
HsCCHCN 
H2CCHCN 
HsCCHCN 
HsCCHCN 
H2CCHCN 
HsCCHCN 
HoCCHCN 
H2CCHCN 
|sCCHCN 
HoCCHCN 
H2CCHCN 
HsCCHCN 
HoCCHCN 
HsCCHCN 
HsCCHCN 
HoCCHCN 
HsCCHCN 
HsCCHCN 
HoCCHCN 
HsCCHCN 


oRoRoRoRoRoReRogoRekogoRekonoRe Rooke koRoke ke) 


H3CgN 
HsCaN 
HsCaN 
H3C3N 
HsCaN 
HsC3N 
HsC3N 
HsCaN 
H3C3N 
HsCaN 
HsCaN 
H3C3N 
HsCaN 
gCgN 
H3C3N 
HsCaN 
H3CgN 
H3C3N 
HsCaN 
H3C3N 
H3C3N 
HsCaN 
H3CgN 
HsCaN 
HaC3N 
HaCaN 


OP PEO EP EEE Eee Eee Eee eee eeeEeEe?e 


Frequency (GHz) 


31.6156000 
35.3790500 
36.0646889 
36.1402200 
36.2225000 
36.8784700 
40.4257132 
41.1870836 
41.2976561 
41.4207138 
42.1384521 
45.4695206 
46.2978836 
46.4528399 
46.6280696 
47.3948488 
75.6463060 
75.6587572 
76.6665519 
78.0996741 
78.8285795 
81.674936 
84.0442448 


31.8489874 
33.8637212 
37.1392132 
38.1027037 
38.3423461 
38.6167091 
39.4865865 
42.0111351 
42.4217854 
43.8024266 
44.2102030 
45.0990808 


33.0513004 
33.0516190 
37.1826557 
37.1830142 
41.3139909 
41.3143891 
45.4453036 
45.4457416 
49.5765916 
49.5770694 
74.3636580 
74.3643850 
82.6257600 
82.6265180 
86.756698 
86.757524 

90.8876208 
90.8884956 
95.0184892 
95.0194037 
99.1493060 
99.1502601 
103.280069 
103.2810626 
107.4107758 
107.4118089 


“The format of the quantum numbers is Jk, 


Transition? 


73,4 — 82,7 
15445 - 14444 
150,15 — 140,14 
159.43 — 149,12 
18444 - 144.43 
160,16 — 150,15 
161,45 — 15444 


for HCCCH2CN and CHsCCHCN (asymmetric rotors), 


logAu (s~*) 


—5,.5633 
—5.4086 
-5.3748 
—5.4089 
—5.4059 
—5.3546 
—5.2292 
—5.1981 
5.2225 
—5.2187 
-5.1751 
—5.0716 
—5.0429 
—5.0603 
—5.0555 
5.0176 
—6.2551 
4.3956 
-4.3767 
—4,3600 
-4.3421 
—4.2935 
—4.2576 


—-5.2773 
-5.1974 
—5.0691 
5.0271 
—5.0555 
—5.0461 
—4,9892 
-5.4001 
—4.8900 
—4,8695 
-4.8574 
—4.8103 


5.3573 
—5,.3505 
-5.1995 
—5.1942 
—5.0590 
—5.0546 
—4,9321 
—4,9286 
-4.8166 
-4.8136 
—4,2809 
4.2795 
—4,1422 
—4.1411 
—4.0780 
4.0770 
—4.0169 
—4.0159 
-3.9584 
—3.9576 
—3.9026 
-3.9017 
—3.8490 
—3,8482 
-3.7975 
—3.7968 


and Jx for CH3C3N (symmetric top molecule). 


Eu (K) 


6.4 
7.9 
6.9 
11.4 
11.4 
8.2 
9.9 
8.9 
13.4 
13.4 
10.2 
12.0 
144 
15.6 
15.7 
12.5 
17.0 
30.2 
29.6 
34.4 
31.4 
33.5 
35.5 


6.2 
6.5 
8.0 
7.3 
10.6 
10.7 
8.4 
4.6 
10.0 
12.7 
12.8 
10.6 


14.6 
7A 
16.4 
8.9 
18.4 
10.9 
20.6 
13.1 
23.0 
15.5 
41.4 
33.9 
49.1 
41.6 
53.3 
45.8 
57.7 
50.2 
62.2 
54.7 
67.0 
59.5 
71.9 
64.4 
774 
69.6 


Blending 


blended with aGg’-(CH2OH)2 
blended with CHzCOCH3 


unb 
unb 
unb 


lended 
lended 
lended 


blended with U 


unb 


lended 
lended 
lended 
lended 


lended wi 
lended wi 
inblended 
lended wi 
lended wi 
inblended 


lended wi 
lended wi 
lended wi 
lended wi 


h c-C3H2 
h CHsNHCHO 


h CHgCONH2 
h HCzCHO 


hU 
hU 
h CoaHsO 
hU 
h CH3COOH 


piss 


unb 


b 
b 
ul 
b 
b 
ul 
blended wi 
b 
b 
b 
b 
b 


lended with C2H<O! 


pin 


lended 


lended 
lended 
lended 
lended 
lended 
lended 
lended 


blended with n-C3H7CN 


unb 


lended 


blended with HNCO 


unb! 
unb 


unb 


lended 
lended 


lended 


blended with CHzO0CHO 
blended with aGg’-(CH2OH)2 
unblended 

unblended 

unblended 

blended with s-C2HsCHO 
blended with H'SCO* 
blended with CHsCOCH3 
unblended 

unblended 

unblended 


blended with CHsCONH» 
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TABLE 4 | Derived physical parameters of the nitriles towards G+0.693 analysed in this work using MADCUBA, along with their associated uncertainties. The fixed 


parameters used in the fit are shown without associated uncertainties (see text). 


Molecule N? (x 10°? cm-?) Tex (K) 
HOCN 2.13+0.04 7,440.2 
HCOCN 0.76+0.11 7A 
HOCH2CN 0.8+0.2 7A 
HCOCH2CN <3.6 7A 
CH2CCHCN 2.3440.06 12.1+0.5 
HCCCH2CN 1.77+0.08 12.1 
CH3CCCN 1.35+0.03 18.6+1.0 


Visr (km s~') FWHM (km s~') Abundance? (x10-'°) 

68,0+0.2 19.2+0.3 16 

68 19.2 0.6 

67 19.2 0.6 

67 19.2 <2.7 
66.1+0.3 21.3+0.7 17 
67.0+0.6 21.3 1.3 

68 21.3 1.0 


“The uncertainties of the column densities are derived by the AUTOFIT algorithm implemented in MADCUBA (see Martin et al., 2019 for details), and do not contain calibrations errors, 


which are expected to be ~ 10%. 
’We adopted N(H2) = 1.35 x 10°? cm”, from Martin et al. (2008). 


Table 4). The rest-frequencies have then been computed in the J = 
0-70 interval with Ky ina, = 50. Theoretical values of dipole 
moments fg = 0.932 D, uy = 1.574 D, and p, = 1.274 D, computed 
at CCSD level (Mollendal et al., 2012) have been employed. All 
the calculations have been performed with the CALPGM suite of 
programs Pickett (1991). 

This species is not detected towards G+0.693. We have derived 
an upper limit for its abundance using the brightest transition 
according to the LTE model that are unblended, namely the 
675-5),4 transition at 101.598576 GHz. MADCUBA calculates 
the upper limit of the column density using the 30 value of the 
integrated intensity (see details in Martin et al., 2019). We have 
used the same T.,, FWHM, and vysp used for HOCH,CN. We 
obtained an upper limit of the HCOCH,CN abundance of 
<2.7 x 107!° (Table 4). 


3.2 C,H3N Isomers 


We report in this section the first detection towards G+0.693 of 
cyanoallene (CH,CCHCN), propargyl cyanide (HCCCH,CN), 
and cyanopropyne (CH3;CCCN). 


3.2.1 Cyanoallene (CH2CCHCN) 


Figure 5 shows the molecular transitions of CHyCCHCN that are 
unblended, or only slightly blended with other species already 
identified in this source, whose spectroscopic information is 
presented in Table 3. The G+0.693 cloud is the second 
interstellar source where CH,CCHCN has been detected, after 
the cold cloud TMC-1 (Lovas et al., 2006; Marcelino et al., 2021). 
We left N, T.., FWHM, and vy sp as free parameters, and obtained 
a column density of (2.34+ 0.06) x 10 cm’, and a molecular 
abundance of 1.7 x 107'° (Table 4). 


3.2.2 Propargyl Cyanide (HCCCH2CN) 


Figure 6 shows the molecular transitions of HCCCH2CN that are 
unblended, or only slightly blended with other species already 
identified in this cloud, whose spectroscopic information is 
presented in Table 3. As in the case of its isomer 
CH,CCHCN, G+0.693 is the second interstellar source where 
HCCCH,CN has been detected, after the cold cloud TMC-1 
(McGuire et al., 2020; Marcelino et al., 2021). We fixed T,,. and 
FWHM to the values obtained for CH,CCHCN, and left N and 


Vispr free. We obtained a column density of (1.77+ 0.08) x 10! 
cm’, and a molecular abundance of 1.3 x 107°°. The 
CH,CCHCN/HCCCH,CN ratio is ~ 1.3. 


3.2.3 Cyanopropyne (CH3CCCN) 


Figure 7 shows the spectra of multiple unblended or slightly 
blended transitions of CH;CCCN (listed in Table 3). Unlike its 
isomers, which are asymmetric molecules, CH;CCCN is a 
symmetric top molecule. For the analysis, we have used the 
lowest energy K = 0 and K = 1 transitions (see Table 3), which 
are the ones that dominate the line emission in a source with low 
Tex like G+0.693 (5 — 20 K; see e.g. Zeng et al., 2018). We fixed the 
FWHM and vygp to the values derived for CH,CCHCN, leaving N 
and T.x as free parameters. We obtained a column density of 
(1.35+ 0.03) x 10!3 cm” (Table 4), and a molecular abundance of 
1.0 x 107!°. The isomeric ratios of CH;CCHCN/CH;CCCN and 
HCCCH,CN/CH3CCCN are ~1.8 and ~1.3, respectively. 


4 DISCUSSION 


4.1 Interstellar Chemistry 
4.1.1 Oxygen-Bearing Nitriles 


We show in Figure 8 the molecular abundances of the O-bearing 
nitriles detected towards G+0.693 studied in this work. The 
relative ratio of the detected species HOCN:HCOCN: 
HOCH,CN is 2.8:1:1. By extrapolating the hydroxy/aldehyde 
(OH/HCO) ratio of HOCN/HCOCN to HOCH,CN/ 
HCOCH,CN, one should expect an abundance of 0.15 x 107'° 
for HCOCH2CN, more than one order of magnitude lower than 
the upper limit derived from current observations (< 2.7 x 1071, 
see Table 4). This suggests that deeper observations reaching 
higher sensitivity will be needed to address the detection of this 
species. 

In the following, we discuss possible formation routes of the 
different O-bearing nitriles, combining the results obtained in 
G+0.693 and in other interstellar sources with theoretical and 
experimental works: 

e HOCN: besides G+0.693, this species was detected 
previously towards several other positions of the Sgr B2 region 
in the Galactic Center (Briinken et al., 2009; Briinken et al., 2010), 
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FIGURE 6 | Selected transitions of propargyl cyanide (HCCCH2CN; see Table 3) detected towards the G+0.693 molecular cloud. The best LTE fit derived with 
MADCUBA for the HCCCH2CN emission is shown with a red curve, while the blue curve shows the total emission considering all the species identified towards this 
molecular cloud. The y-axis shows the line intensity in antenna temperature scale (7,4) in Kelvin, and the x-axis shows the frequency in GHz. 
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FIGURE 7 | Selected transitions of cyanopropyne (CH3sCCCN; see Table 3) detected towards the G+0.693 molecular cloud. The best LTE fit derived with 
MADCUBA for the CH3CCCN emission is shown with a red curve, while the blue curve shows the total emission considering all the species identified towards this 
molecular cloud. The y-axis shows the line intensity in antenna temperature scale (7,4) in Kelvin, and the x-axis shows the frequency in GHz. 
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HOCN 


FIGURE 8 | Molecular abundances with respect to Hs of the oxygen-bearing nitriles studied in this work derived in different interstellar sources. Purple bars 
correspond to G+0.693 (this work; see Table 4), with the HCOCH2CN value indicating an upper limit. We compare with other sources: several positions also in the Sgr 
B2 region (magenta; Brunken et al., 2010, see also Marcelino et al., 2010); several dense cores (B1-b, L1544, L183, and L483) and the lukewarm corino L1527 (green: 
Marcelino et al., 2010; Marcelino et al., 2018); the dark cloud TMC-1 (yellow; Cernicharo et al., 2021); and the IRAS 16293-2422 B hot corino (red, Zeng et al., 
2019). To derive the uncertainties of the molecular abundances we have considered the uncertainties of the molecular column densities reported in the different works, or 
a 15% of the value of N if the uncertainty is not provided, and we assumed an uncertainty 
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FIGURE 9 | Upper panel: Abundances with respect to He of the CsH,N isomers detected towards G+0.698 (purple; this work) and TMC-1 (yellow, Marcelino et al., 
2021). To derive the uncertainties of the molecular abundances we have considered the uncertainties of the molecular column densities of the C3H,4N isomers reported in 
this work (Table 4) and in Marcelino et al. 2021, and we assumed an uncertainty for the (Hz) column density of 15%. Lower panel: Molecular ratios between the 


CH,CCCN 


and towards several dense cores (B1-b, L1544, L183, and L483) as 
well as the lukewarm corino L1527 (Marcelino et al., 2010; 
Marcelino et al, 2018). Figure 8 shows that the HOCN 
abundance derived in G+0.693 is of the same order of 
magnitude of those detected in other Sgr B2 positions 
(~ 10°" — 1071; Briinken et al., 2010), and higher than those 
derived in the dense cores and L1527 (Marcelino et al., 2018). 
This suggests that the role of surface-chemistry and the presence 
of shocks enhance the HOCN abundance, similarly to its 
isomer HNCO (Hasegawa and Herbst 1993; Garrod et al., 
2008; Martin et al., 2008; Rodriguez-Fernandez et al., 2010; 
Quénard et al., 2018). The chemistry of the molecular clouds of 
the Galactic Center, and that of G+0.693 in particular, is 
dominated by large-scale shocks (Martin-Pintado et al. 


2001; Martin et al., 2008), which are responsible for the 
sputtering of dust grains, releasing many molecules formed 
on the grain surfaces into the gas phase (see Caselli et al., 1997; 
Jiménez-Serra et al., 2008). This can increase the abundance of 
the species by orders of magnitude. Similarly to isomer HNCO, 
which is efficiently formed on grain surfaces by hydrogenation 
of accreted OCN (Hasegawa and Herbst 1993; Garrod et al., 
2008), HOCN can also be formed on grain mantles if the 
oxygen atom is hydrogenated: 


OCN + H > HOCN, (1) 


and then subsequently released by shocks (Briinken et al., 
2010). An alternative surface route might be the reaction of two 
highly abundant species: 


Frontiers in Astronomy and Space Sciences | www.frontiersin.org 


13 


July 2022 | Volume 9 | Article 876870 


Rivilla et al. 


New Nitriles in G+0.693-0.027 


CN + OH — HOCN (2) 


e HCOCN: this species has been detected previously in the 
massive hot core SgrB2 (N) (Remijan et al., 2008), and in the dark 
cloud TMC-1 (Cernicharo et al., 2021). The HCOCN abundances 
found in G+0.693 and TMC-1 are very similar, in the range of 
(3.5-5) x 107!!, as shown in Figure 8. These two regions have 
very different physical conditions, which imprint their chemistry. 
While in the case of the dark and cold TMC-1 cloud gas-phase 
chemistry is thought to be dominant, since thermal or shock- 
induced desorptions are highly unlikely, the chemistry of 
G+0.693 is strongly affected by shocks, and thus surface 
chemistry also plays an important role. Therefore, the similar 
HCOCN abundances in G+0.693 and TMC-1 points towards a 
predominant gas-phase chemistry origin. Indeed, the quantum 
chemical calculations by Tonolo et al. (2020) have shown that 
HCOUCN species can be efficiently formed through the gas-phase 
reaction between formaldehyde (H,CO) and the cyanide radical 
(CN), which are highly abundant species in the ISM, in which the 
CN radical attacks the unsaturated carbon of H,CO and 
substitutes one of the H atoms: 


H,CO + CN — HCOCN + H, (3) 


e HOCH,CN: this species was first detected in the ISM towards 
the hot corino IRAS 16293-2422 B (Zeng et al., 2019), and more 
recently towards the SMM1 hot corino in Serpens (Ligterink et al., 
2021). The abundance derived in G+0.693 is 4.3 x 10°1°, very similar 
to that derived in the hot component of IRAS 16293-2422 B 
(Figure 8). The chemical model by Zeng et al. (2019) considered 
the surface formation route proposed by the laboratory experiments 
of Danger et al. (2012); Danger et al. (2013): 


H,CO + HCN — HOCH,CN (4) 


and ion-neutral destruction reactions with H}, HCO", and H;0°, 
and concluded that more chemical pathways are needed to explain 
the abundance observed in the hot corino IRAS 16293-2422 B. More 
recently, the quantum chemical cluster calculations performed by 
Woon (2021) have proposed new surface reactions between C*, 
which is distributed throughout the whole Galactic Center (Harris 
et al., 2021), and two very abundant species, HCN and HNC (e.g. 
Colzi et al. 2022), embedded in H,O icy grain mantles. The C* ion 
reacts with HCN and HNC forming intermediate species that attacks 
neighboring H,O molecules of the ices, resulting into the radicals 
HOCHNC and HOCHCN. These species can be easily hydrogenated 
on the grain surfaces to form HOCH,CN. The inclusion of these 
alternative surface routes in the chemical models might help to 
explain the HOCH2CN abundances detected in G+0.693 and hot 
corinos, where the molecules can be injected to the gas phase through 
shocks and thermal effects, respectively. 

¢ HCOCH,CN: the theoretical calculations performed by 
Horn et al. (2008) proposed that this species might be formed 


from two abundant precursors in the ISM: 
HC3N + H,O — HCOCH,CN (5) 


However, while this reaction might occur in aqueous solution, 
its activation energy, 216 kJ mol ? (25,980 K), is too high to occur 


in the ISM. Recently, Alessandrini and Melosso (2021) have 
studied the reaction between oxirane (or ethylene oxide, c — 
C,H,O) — also detected towards G+0.693 (Requena-Torres et al., 
2008) — and the CN radical. Although the main pathway is the H 
abstraction from oxirane, forming the oxiranyl radical, the 
formation of HCOCH,CN + H is also possible with a rate of 
~ 10° cm? molec’ s~'. New theoretical and/or experimental 
works of this species are needed to determine if it can be 
efficiently formed in the ISM, opening the possibility for its 
interstellar detection. 


4.1.2 C,H3N Isomers 


The unsaturated C,H3N isomers towards G+0.693 have very 
similar abundances within a factor of 2, spanning a range of 
(1.0-1.7) x 10°!° (Table 4), as also previously observed in the 
dark cloud TMC-1 by Marcelino et al. (2021). Moreover, Figure 9 
shows that the abundances in these two molecular clouds, which 
have very different physical conditions, as mentioned above, are 
very similar. This suggests that these molecules are 
predominantly formed through gas-phase chemistry (see 
previous discussion about HCOCN). Furthermore, since the 
three isomers are almost equally abundant, their respective 
formation might be linked to common precursors. Indeed, 
Balucani et al. (2000) proposed that these unsaturated nitriles 
can be formed efficiently by reactions in which the cyanide radical 
(CN) attacks an unsaturated carbon of the hydrocarbons 
methylacetylene (CH3;CCH) and allene (CH,CCH,): 


CH;CCH + CN 222222, CH,CCCN + H, (6) 
ah CH,CCHCN +H, (7) 

CH,CCH, + CN 2°8CH,CCHCN + H, (8) 
°*, HCCCH,CN + H; (9) 


in which the branching ratios for each reaction are indicated 
above each arrow (normalized to 1). These ratios were derived 
using the experiments and quantum chemical calculations by 
Abeysekera et al. (2015)/Balucani et al. (2000) in the first two 
reactions, and from Balucani et al. (2002) in the latter two 
reactions. These radical-neutral reactions show no entrance 
barriers, they have exit barriers well below the energy of the 
reactant molecules, and are exothermic. The proposed precursors 
CN and CH3CCH are abundant molecules in the ISM. In 
particular, they were detected towards G+0.693 with molecular 
abundances of 1.5 x 10° and 1.3 x 10°%, respectively (Rivilla et al., 
2019a; Bizzocchi et al., 2020), so they are viable precursors. Allene 
(CH,CCH2) has zero dipole moment, so its detection through 
rotational spectroscopy is not possible, and thus its abundance is 
unknown. However, the similar abundances of the three isomers 
suggest that it can be as abundant as CH3CCH in the ISM. 

Regardless of the actual abundance of CH,CCH3, which is 
unknown, the proposed branching ratios seem to be in conflict 
with the observational findings in G+0.693 and TMC-1, since 
they are not able to produce equal abundance for the three 
isomers. As already noted by Marcelino et al. (2021), it would 
be interesting to study the branching ratios of the CH,CHCH) + 
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CN reaction using the chirped-pulse uniform flow experiment 
used by Abeysekera et al. (2015) for the CH,CCH + CN reaction, 
and compare them with the values derived from quantum 
chemical calculations by Balucani et al. (2002), to reconcile the 
experimental/theoretical works with the findings of the 
observations. 


5 CONCLUSION: IMPLICATIONS FOR THE 
RNA-WORLD 


Compounds of the nitrile family, under early Earth conditions, offer 
a rich chemistry due to the large number of reactions that they can 
trigger. Nitriles could be transformed into amides, carboxylic acids 
and esters via hydrolysis and alcoholysis _ respectively. 
Autocondensation of nitriles in a basic environment could yield 
to cyanoketones and cyanoenamines, a high reactive intermediate in 
the synthesis of complex five- and six-member heterocycles (Erian 
1993). The high amounts of ammonia of the reducing atmosphere of 
the primitive Earth is a favorable scenario to obtain amidines from 
nitriles (Shriner and Neumann 1944). Moreover, the NCN backbone 
of amidines offer an unique structure to yield complex N-containing 
heterocycles like purine and pyrimidine nucleobases. Furthermore, 
nitriles can activate the formation of the building blocks of RNA, 
ribonucleotides (e.g. Powner et al., 2009; Powner and Sutherland 
2010; Patel et al., 2015). Two of the nitriles studied in this work, ie. 
glycolonitrile and cyanoacetaldehyde, have been proposed as 
activation agents for the formation of more complex molecules 
with prebiotic relevance. The latter (HCOCH,CN) is a precursor of 
cytosine (Robertson and Miller 1995; Nelson et al., 2001; Menor- 
Salvan et al. 2009). The former (HOCH,CN) is not only a 
fundamental precursor to ribonucleotides and lipids (Ritson and 
Sutherland 2012, 2013; Patel et al. 2015; Liu et al. 2018; Ritson et al. 
2018), but also of other biologically-important molecules such as the 
simplest amino acid glycine (NH,CH,COOH; Rodriguez et al., 
2019), and of the nucleobase adenine through rapid HCN 
oligomerisation (Schwartz and Goverde 1982; Menor-Salvan and 
Marin-Yaseli 2012). Unsaturated carbon-chain nitriles like the 
C4H3N isomers studied in this work are also especially interesting 
for prebiotic chemistry because the presence of unsaturated bonds 
allows further chemical evolution that can produce biomolecules 
(Rosi et al. 2018). 

This work extends the repertoire of nitriles detected in the 
G+0.693 molecular cloud, a region that exhibits one of the richest 
chemical content in the ISM, and hence it is a well suited testbed 
to census the molecular species present in the ISM. Besides 
HOCN, already reported by Briinken et al. (2010) and Zeng 
et al. (2018), we have provided the tentative detections towards 
this source of HCOCN and HOCH,CN (third detection in the 
ISM), and the detection of the three unsaturated C,H3N isomers 
(being the second source after TMC-1 in which all three isomers 
are identified). These detections confirm the rich reservoir of 
nitriles in space, and complete the list of prebiotic molecular 
precursors detected previously, including species directly 
involved in the synthesis of ribonucleotides such as 
glycolaldehyde (HCOCH,OH; Hollis et al. 2004; Requena- 
Torres et al., 2006; Beltran et al., 2009; Jorgensen et al., 2012), 


urea (Belloche et al., 2019; Jiménez-Serra et al., 2020), 
hydroxylamine NH,OH (Rivilla et al. 2020), and 1,2- 
ethenediol (Rivilla et al. 2022a); of amino acids, such as amino 
acetonitrile (NH,CH,CN; Belloche et al., 2008; Melosso et al., 
2020); and of lipids, such as ethanolamine (NH,CH,CH,OH; 
Rivilla et al., 2021a), and propanol (CH;CH,CH,OH; Jimenez- 
Serra et al. 2022; Belloche et al. 2022). 

In star- and planet-forming regions, this chemical feedstock 
can be processed through circumstellar disks, and 
subsequently incorporated into planetesimals and objects 
like comets and asteroids. We know that our planet suffered 
a heavy bombardment of extraterrestrial bodies ~500 Myr after 
its formation (e.g. Marchi et al., 2014). Laboratory impact 
experiments have shown that a significant fraction of the 
molecules contained in comets and meteorites could have 
been delivered intact to the early Earth (Pierazzo and 
Chyba 1999; Bertrand et al., 2009; McCaffrey et al., 2014; 
Todd and Oberg 2020; Zellner et al., 2020). Once on the 
planetary surface, under the appropriate physical/chemical 
conditions, these molecules could have allowed the 
development of the prebiotic processes that led to the dawn 
of life on Earth. 
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